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Magnetic Sd-ions doped into wide-gap oxides show signatures of room temperature ferromag- 
netism, although their concentration is two orders of magnitude smaller than that in conventional 
magnets. The prototype of these exceptional materials is Co-doped ZnO, for which an explana- 
tion of the room temperature ferromagnetism is still elusive. Here we demonstrate that magnetism 
originates from Co^^ oxygen-vacancy pairs with a partially filled level close to the ZnO conduction 
band minimum. The magnetic interaction between these pairs is sufficiently long-ranged to cause 
percolation at moderate concentrations. However, magnetically correlated clusters large enough to 
show hysteresis at room temperature already form below the percolation threshold and explain the 
current experimental findings. Our work demonstrates that the magnetism in ZnOiCo is entirely 
governed by intrinsic defects and a phase diagram is presented. This suggests a recipe for tailoring 
the magnetic properties of spintronics materials by controlling their intrinsic defects. 



ZnO is a piezoelectric conductive oxide, in which free- 
carriers coexist with optical transparenc}^^ ^. If made 
magnetic, ZnO will become the ultimate multifunctional 
material, with semiconducting, magnetic, optical and 
mechanical properties. This will have a far reaching 
impact on the emerging field of spintronics^ with ap- 
plications in optoelectronic J2 and quantum computing^ . 
Moreover it will allow us to go beyond the (Ga,Mn)As 
paradigm^, whose practical use is severely hampered 
by the low ferromagnetic critical temperature. This is 
why ZnO: Co is perhaps the most studied among all the 
diluted magnetic oxides. Room-temperature ferromag- 
netism (RTF), first demonstrate d hj U eda et al.^, is now 
confirmed by a number of groups^ ^ (see Table I of the 
supplementary material). The experimental situation is 
however still confused and here we list the main findings. 

i) Spectroscopy confirms that Co^+ substituting Zn 
is the center responsible for all the differe nt mag- 
netic phases found experimentally, including RTF*^^^^^^, 
paramagnetism^^^^ and spin-glasseJ^^. RTF is usually 
assigned from magnet ometr};^ ^ ^ 

ii) The saturation magnetization Mg and the rema- 
nence are always small and secondary phases are of- 
ten difficult to rule out. However, except for metallic 
Co, most of them are either non-magnetic or antifer- 
romagnetic with low Neel temperatures (CoO, C02O3, 
C03O4, ZnCo204). The coercive field is typically small 
(^100 Oe) and only weakly temperature dependent. 

iii) Mg is usually smaller than what is expected 
for Co^+ with values as low as 0.01 /ie/Co^^, sug- 
gesting antiferromagnetic interaction among Co^+ and 
frustration^^ 

iv) Growth conditions and annealin g are crucial for 
the magnetic state. Chemical methods 1 ^^ * and molec- 
ular beam epitax};^^ generally lead to paramag netism, 
while pulsed laser deposition produces RTF film d^ l ^ * ^ ^. 
Typically oxygen deficient growth^^ at tuned substrate 
temperatures^^ promotes RTF. Similarly, annealing 
in vacuum enhances the magnetic moment and produces 
ferromagnetism^^ while annealing in oxygen has 



the opposite effeclR*^. 

v) The role played by free carriers in establishing mag- 
netism is unclear. Sequential annealing in reducing and 
oxidizing atmosphere reveals little correlation between 
the electrical conductivity and the magnetic statd^. 
Similar conclusions are reached for Al- and H-doping'^. 

vi) ZnO is often reported to be oxygen deficient. This 
has been attributed to either oxygen vacancies 

or multi-center H^^, with Zn interstitials ( Zn^) n ow ruled 
out by both experimental^^ and theoreticaP^I^ evidence. 
Thus the promotion of RTF due to Zn vapour exposur^^ 
cannot be attributed to an increase of the Zn^ concentra- 
tion. 

vii) Electron paramagnetic resonance^^ suggests the 
presence of two magnetic centers. These are both re- 
lated to Co^+, although they exhibit fine differences in 
the signal. Interestingly, for a nominal Co concentration 
about 5% the two centers appear with similar abundance. 

Existing mechanisms for ferromagnetism in the diluted 
limit cannot explain this complex collection of phenom- 
ena. The p-d Zener model^^ lacks of its foundations when 
applied to Zni_a;Coa;0 since t here is little correlation be- 
tween carriers and magnetisnP^l^. When present, carri- 
ers are electrons and not holes with small exchange cou- 
pling to the local spins and therefore the typical critical 
temperatures (Tc) are tiny at realistic carrier densities. 
Similarly super-exchange must be ruled oulP^. This is 
short ranged and RTF can be obtained only for x above 
the nearest neighbour (NN) percolation threshold. For 
the wurtzite lattice this is 20 %, much greater than the 
typical experimental concentrations. 

Finally a modification of the Zener scheme, called the 
"donor impurity band exchange" (DIBE), assumes that 
the magnetic interaction is mediated by large hydrogenic 
orbitals associated to intrinsic defects and predicts fer- 
romagnetism below the donor percolation threshold, i.e. 
in absence of free carriers^^. Although the mean field Tq 
obtained with realistic parameters for Zni-^^Coa^O is ex- 
tremely small (~10 K)^, the model is frequently used to 
explain the experimental results. We have investigated 
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such a model with Monte Carlo calculations (see supple- 
mentary materials) and demonstrated that for Co and 
donor concentrations respectively of 10 % and 1 %, the 
Tc is only a few degrees K. RTF is obtained only with 
unrealistically high values of the exchange coupling, and 
therefore the DIBE scheme must also be ruled out. 

In the absence of any simple-scheme for ferromag- 
netism we turn to atomistic density functional theory 
(DFT) . For this problem the standard approximations to 
the exchange and correlation potential (LDA and GGA) 
are not appropriate since they over-delocalize and under- 
bind the Co-d shell resulting in their incorrect positioning 
with respect to the Fermi level, E^. These failures are 
only minor in the case of GaAsiMiP^, but they become 
a serious drawback for the oxides where the Co d density 
of states (DOS) has important contributions in the ZnO 
band-gap. We therefore use an approximated version of 
the self-interaction correction (ASIC) scheme^^, which is 
free from these problems and produces exchange param- 
eters for transition metal monoxides superior to those 
obtained with LDA/GGA^. 

In Fig. fl] we present the DOS of a Co impurity (at 
the Zn site) in a 128 atom ZnO supercell (x=0.0156) cal- 
culated with both LDA and ASIC. Although they both 
predict a 2+ valence, the position of the Co d levels is 
remarkably different in the two cases. LDA places the 




FIG. 1: Density of states of Zni-xCoccO as calculated from 
density functional theory. The simulation cell is a 128 ZnO 
supercell in which one Zn atom is replaced by Co (x=0.0156). 
Panels (a) and (b) show LDA and ASIC results respectively. 

occupied minority e states just below at the edge of 
the ZnO conduction band minimum (CBM). ASIC shifts 
these by about 2 eV down to the valence band top (VBT) 
in agreement with recent calculations^^. In addition the 
ZnO band-gap opens and the Zn d DOS is also down- 
shifted. This gives us a picture where there is no Co-d 
contribution to the DOS around E-p^ with the first unoc- 
cupied minority states (^2) placed at about 1 eV above 
the CBM. Such a DOS is in excellent agreement with 
ultraviolet photoemission (UPS), which places the Co d 
DOS at a binding energy of about 3 eV with a satellite 
peak at 7 eV and a rather diffuse taiP^. This electronic 



structure is incompatible with a carrier mediated picture 
of ferromagnetism. In fact, E^ should be moved by at 
least 1 eV in order to affect the Co valence and there- 
fore to promote the charge transfer necessary for strong 
exchange coupling^^, a task hardly achievable. 

We test this conjecture by calculating the magnetic 
coupling between two Co^+ at various distances dco- Co- 
in table |l] we show the magnetic coupling energy Em 
{Em > indicates ferromagnetic coupling), obtained 
as the total energy difference between the ferromagnetic 
(FM) and antiferromagnetic (AFM) configuration of a 
selection of supercells containing two Co^+ and one in- 
trinsic defect far from both the Co^+. Clearly the ex- 
change interaction between two Co^+ is strong only at 
NN position, i.e. when the super-exchange interaction 
is effective. In this case the interaction is AFM in 
the a-b plane (dco-Co=3.19A) and FM along the c axis 
(dco-Co=3.11A). Em however drops to zero already at 
second NN, regardless of the presence of additional in- 
trinsic defects. In particular we show data for third NN 
(the data for second NN are similar) in presence of either 
Zn^, H or Vq from which one has to conclude that RTF 
is not achievable by simply defect doping. 

Figure [2] offers an insight on why Zn^ and Vq are un- 
able to mediate RTF. We present the DOS of a 128 atom 
supercell containing one Co^+ and one intrinsic defect 
(Vo, Zn^ and H), and compare the case where the Co^+ 
and the defect are well separated in the cell with that 
in which they are at NN position. In the case of distant 
defects the DOS is essentially a superposition of that of 
Co^+ and the defect. Both Zn^ and H posses a filled hy- 
drogenic level above the CBM and rather close to the Co 
d-t2 minority levels. In contrast Vq displays a doubly oc- 
cupied impurity level 1 eV above the VBT, almost at the 
same position of the minority Co-d e levels. Most impor- 
tantly there is no evidence of interaction between Co^+ 
and the defect levels (the magnetic moment calculated 
from the Mulliken population is ~2.6 /in/Co for all three 
cases, similarly to the case of Co^+ only). This means 
that, despite the energy proximity of the defect levels to 
those of the Co 3d shell, the overlap of the hydrogenic 
wave-function at the Co site is only minimal. For this 
reason we conclude that the DIBE model as formulated 
cannot be sustained by the electronic structure of Co^+ 
in ZnO. 

The right panels of Fig. |2] for NN defects give a dif- 
ferent picture. In the case of both H and Zn^ there is a 
substantial charge transfer from the defects to the Co^+ 
resulting in a partial occupation of the minority ^2 levels 
and a reduction of the magnetic moment (1.95 ji^ /Co 
and 2.05 (jlb/ Co respectively for Zn^ and H). On-site re- 
pulsion moves upwards in energy all the Co 3d manifold 
and now the majority ^2 levels occupy the ZnO bandgap. 
In contrast, for Vq slight charge transfer from the Co 
3d to the vacancy and hybridization move the 3d man- 
ifold downwards in energy and the doubly occupied Vq 
level spin-splits in the opposite direction to that of Co. 
Such a feature suggests that the exchange interaction be- 
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TABLE I: Calculated magnetic energy Em for various magnetic centers and different dopants configuration. CI and C2 are 
the two magnetic centers included in the simulation cell and their relative concentration (concentration of each center), D is 
the dopant with its concentration, dci-C2 in the distance between the two center expressed both in A and in NN shells. For 
some of the NN complexes we present the geometrical configuration (after relaxation in the pictures below). 
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tween NN Co ions, when mediated by a defect, can be 
extremely large. This is indeed the case as demonstrated 
in table [l] where Em for various Co-Co-defect complexes 
is presented. Such complexes effectively behave as small 
metallic clusters and it is of no surprise that the exchange 
energy increases quite dramatically for dco-Co around 
and below 2.5 A, which is the NN distance in metallic 
Co. 

Is this sufficient for RTF in the diluted phase? Un- 
fortunately not. These interactions, although strong, are 
still short-range and therefore produce RTF only for x 
around the percolation limit. Moreover, the strong FM 
interactions require a ratio between Co and the relevant 
donor of 2:1, which at percolation means typical defect 
concentrations of around 10%. These are impossible to 
achieve for any of the defects investigated even under 
the most favorable conditions. Finally those complex 
structures, if abundant, should appear spectroscopically 



with a substantial increase of the Co-d DOS in the ZnO 
bandgap. 

Since Co^+ alone cannot be responsible for RTF at 
low dilution we have searched for other possible magnetic 
centers and found that Co^+-Vo pairs (CoV) overcome 
the limitations mentioned above. In figures |2]i and |2^ 
we present the DOS associated to the Co-d shell for both 
Co^+ and CoV, as compared with UPS from referencd^. 
Clearly both Co^+ and CoV are compatible with UPS, 
in particular they both show a finite DOS at about -8 eV 
from Ep. This feature is absent in the DOS of both the 
Co-Zn^ and Co-H complexes, which instead present sub- 
stantial contributions in the ZnO gap and therefore are 
incompatible with the spectroscopy. Moreover CoV is 
the only complex among the ones studied which main- 
tains Co in the 2+ valence state. 

CoV are also likely to be abundant. Our calculated Vq 
formation energy (see supplementary material) in zinc 
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FIG. 2: Density of states for a ZnO 128 atom supercell with 
one Co^^ and one additional defect: a) Zn^, b) Vo, and c) H. 
The left panels are for Co^^ and the defect well separated in 
the cell, while the right panels are for the NN position. The 
arrows indicate the relevant defect position. In d) and e) we 
present the Co 3 d density of states as compared with UPS 
data from reference^^. d) Co^^ substitutional at the Zn site, 
e) Co^^-Vo complex. The UPS signal has been aligned to the 
calculated DOS in order to have the first peak at the minority 
Co e states. 



rich conditions is 0.65 eV, which suggests that the Vo 
concentration in ZnO can be as large as 1 % at equilib- 
rium. Since Co does not introduce doping its presence 
will not change considerably the Vo formation energy. 
We have then only to establish whether the Vo sits pref- 
erentially close to a Co site. By using a 128 atom su- 
percell we calculate a reduction in total energy of about 
340 meV for Co and Vo moving from third to first NN 
(pairing energy). This is quite a large gain suggesting 
that most of the Vo are indeed likely to be close to Co 
ions. The large pairing energy also means that the rela- 
tive concentration of CoV with respect to that of 
Co^+ will increase upon oxygen absorbing process- 
ing. This is the case for both annealing in an oxygen poor 
atmospheres and long exposure to Zn and Ti vapours'^. 
In the first case one expects Vo migration to the Co sites, 
while in the second a preferential Vo formation close to 
Co. Finally we found that the pairing energy between 
two Co^+ is also large (~ 210 meV), while we do not find 
any substantial pairing interaction among the CoV in 
any charging configuration. We then expect an inhomo- 
geneous distribution of the Co^+, which in turn would 
lead to an inhomegeneous distribution of CoV and the 
formation of both high and low concentration regions. 



Finally we have to establish whether CoV couple at 
long-range. Table [l] shows Eu for various CoV combina- 
tions. Similarly to the case of Co^+ also the magnetic 
coupling between two charge neutral CoV is remarkably 
weak already at second NN (-6 meV). However, when 
an additional electron donor is present, the situation 
changes dramatically with second and third NN mag- 
netic coupling reaching up respectively to Em ~80 meV 
and Em ~20 meV. Note that we cannot exclude an even 
longer interaction range, which however is hardly acces- 
sible in our simulations since the size of the supercells 
becomes prohibitively large. 
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FIG. 3: CoV impurity band, (a) level diagram for Co^^, Vo 
and CoV, (b) Em for two 2^^ NN CoV as a function of the 
donor impurity band electron density. The electron density 
6.7-10^° cm ~^ corresponds to one electron every two CoV. 

Why does CoV sustain long range coupling while Co^+ 
does not? The cartoon of Fig. [3^ shows the electronic 
structure of Co^+, Vo and CoV. The main feature is 
that the Co^+ empty t2 minority state broadens in CoV 
and forms a t2-Vo hybrid level right at the CBM. This 
extends in space over the the first shell of Zn ions and an 
impurity band forms already at tiny concentrations. Such 
a band can be easily n-doped. Thus, while for Co^+ the 
only exchange mechanisms available either involve virtual 
transitions or weak s-d'^^ exchange interaction with hy- 
drogenic orbitals, for CoV strong carrier-mediated mag- 
netic interaction is possible. In figure ^ we show 
for two CoV in a 128 atom unit cell as a function of the 
electron doping. This is introduced by either moving Ep 
in our simulations or by explicitly introducing an H ion 



in the supercell. Clearly Em depends strongly on the oc- 
cupation of the impurity band, with a maximum at half 
filling (n=6.7xl0^^ cm~^) and vanishing for n=0 (empty 
band). 

Thus DFT offers us a mechanism for the ferromag- 
netism in (Zn,Co)0 based on two magnetic centers. Co^+ 
is responsible only for short range coupling, while CoV 
can instead sustain long range interaction via a fraction- 
ally filled impurity band. Can this alone produce RTF? 
Percolation theorj/^^ sets a strict condition for a magnetic 
ground state of diluted systems: the concentration of 
magnetic impurity should exceed the percolation thresh- 
old Xc. This depends on the range of the interaction and 
for the fee lattice we find 19.8%, 13.7% and 6.2% for inter- 
action extending respectively to 1^*, 2"^ and 3^^ NN^^. 
Therefore our two-center model already produces long 



range ferromagnetism at x 



CoV 



6%. However, we 



do not need such a large x for observing hysteresis in the 
M-H curve at room temperature (i.e. for explaining the 
various experimental claims). This can be achieved below 
Xc since one just needs a number of percolating clusters 
large enough to be superparamagnetically blocked. Note 
that "clusters" here mean regions where x is locally larger 
than Xc The presence of CoV pushes this limit far below 
the 20% needed by Co^+ and by the recently proposed 
model where the magnetism originates from uncompen- 
sated spins at the surface of antiferromagnetic clusters^^. 

The size of those clusters can be estimated by con- 
sidering coherent rotation of the magnetization over an 
anisotropy barrier DNbS^ {D is the zero-field splitting, 
A^B the number of magnetic ions magnetically blocked 
and S the Co spin) . By taking D = 2.76 cm~^ from 
EPR measurements^ we obtaine the estimate A^b ^ 800 
for a blocking temperature Tb = 300 K. This however is 
rather conservative. In granular magnets random dipo- 
lar interaction^^, random magnetic anisotrop}/ or spin- 
odal decompositioiJ^ can push Tb to values considerably 
larger than those predicted for single particle coherent ro- 
tation (up to a factor 5). Thus we estimate A^b in the 
range of 250 magnetic ions. Therefore, an observation of 
a hysteresis requires the existence of regions where 250 
CoV ions interacting at 3^^ NN exist at concentrations 
larger than ~ 6%. Indeed this is a rather modest require- 
ment. 

In Fig. [4^ we present typical cluster distributions P{N) 
as a function of the cluster size N for various concentra- 
tions of Co^+ and CoV. These have been obtained by 
filling randomly a wurtzite lattice comprising 10^ sites. 
As expected P{N) moves from small to large clusters as 



^CoV 



is increased with respect to x . In particular one 
notes that already for = 2% large clusters appear in 
the distribution, which becomes bi- modal at > 6%, 

i.e. above x^ We emphasize that these P{N) have been 
obtained from a completely random distribution, i.e. ne- 
glecting the tendency to clustering suggested by the pair- 
ing energy. 

Finally we investigate the thermodynamical properties 
of our two center model by performing Monte Carlo sim- 
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FIG. 4: Monte Carlo simulations for the two center model: 
(a) cluster distribution P(A^) as a function of the number N 
of magnetic ions in the cluster for different x^^ and x^^^ 
(x^^ /x^^^); (b) magnetization curves m{T) at different Co 
concentration for x^^/x^^^ = 5; (c) m(T) for a total Co 
concentration of 6% and different (d) m(T) for 

x^""^ above percolation (7%) and various (x^°/x^°^). 



ulations for a Heisenberg energy (|S^| = 1) 

Co, CoV Co, CoV 

Heff= J^J^^^^J'^ ^ ^(^^ ' • (1) 



The exchange parameters are chosen to mimic the short 
and long range exchange between Co^+ and CoV respec- 
tively. At NN Jij is AFM for Co'^+ pairs (15 meV) and 
FM for CoV pairs (50 meV) and between Co^+ and CoV 
(50 meV). Moreover it extends to 2"^ (15 meV) and 3'^ 
NN (5 meV) for CoV pairs. The last term accounts for 
an hard- axis easy-plane anisotropy (|n| = ij^^. Note that 
we implicitly assume doping in the CoV impurity band 
and neglect the NN FM component of the exchange be- 
tween Co^+ pairs. Given the uncertainty over the pre- 
cise microscopic configuration and the relative abundance 
of the various complexes (Table |l| our numerical values 
are only representative and certainly conservative. How- 
ever, even with this choice Tq = 250 K above percolation 
j^^CoV _ suggesting that RTF is indeed possible. 

In Fig. (4)3 we present the reduced magnetization m = 
M/Ms as a function of T for different Co concentrations, 
while keeping x^^ /x^^^ = 5. The magnetization curves 
show a transition from a concave upwards shape at low 
concentrations to a convex one for high. As the con- 
centration increases one encounters the two percolation 
thresholds, respectively for CoV and for Co. This pro- 
duces the change in shape, which however is complete 



6 



only above 20%, i.e. when the Co^+'s start to percolate. 
Interestingly m never reaches 1 because of the strong NN 
AFM interaction among Co. 

We also investigate the interplay between Co^+ and 
CoV. In figure Fig. ^ m{T) is calculated by keeping the 
total Co concentration to 6% and by varing the ratio 
between and while in Fig.Jiji we keep x^^^ = 

7% (above percolation) and change x^. In both pictures 
is well below the percolation threshold for NN and 
the high-temperature region of m(T) is almost entirely 
dominated by For instance one may note that for 

^CoV _ 777, (T) approaches zero roughly at the same T 
regardless of whether is 1% (Fig.j^) or 25% (Fig.|4]D 
for a total concentration of 30%). In contrast the amount 
of Co^+ affects the low temperature region, where the 
strong AFM interaction can drastically alter the total 
magnetic moment. For instance m(0) drops from 1 to 
about 0.8 if is increased from 1 to 3%, by keeping 

We are now in a position to propose a phase diagram 
for (Zn,Co)0 based on the relative concentration of Co^+ 
and CoV (Fig. [s]). The most important feature is the 
presence of what we called a blocked superparamagnetic 
phase. This is below and but nevertheless al- 
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probes are needed. In particular further insights would 
be provided by a thorough analysis of small angle neu- 
tron scattering data similar to the case of disordered 
ferromagnets^^. Alternatives may be muon rotation, 
high resolution EPR and energy-dispersive X-ray spec- 
troscopy. 

Finally we partition the long-range FM region into 
two regions separated by the CoV percolation thresh- 
old For > percolation among CoV 
is achieved and one expects measurable conductivity 
from the impurity band. Since the exchange is strong 
an anomalous Hall effect (AHE) should be detected. 
This is not expected for ferromagnetism below 
since the conductivity is then dominated by band con- 
ductivity which is weakly affected by Co^+, given the 
small exchange. Note that this phase diagram says lit- 
tle about the overall conductivity of (Zn,Co)0, which 
in turn can be determined by electrons in the conduc- 
tion band. Moreover, the precise location of the phase 
boundaries depends on details such as the concentration 
of electron donors. In the extreme case of fully compen- 
sated samples the blocked superparamagnetic phase may 
even disappear entirely. 

In conclusion, by using a combination of DFT and 
Monte Carlo simulations, we have demonstrated that 
the observed RTF in (Zn,Co)0 can be attributed to 
blocked superparamagnetism. This develops at concen- 
trations below those required by long-range ferromag- 
netism. However, even this model requires a second mag- 
netic dopant in addition to Co^+ substituting Zn, capable 
of mediating magnetic interaction beyond nearest neigh- 
bours. We have identified the Co-Vq pair as the most 
likely candidate and demonstrated that such center can 
indeed promote long range coupling, if additional n dop- 
ing is present. These findings draw a new roadmap for 
designing diluted magnetic oxides, where the engineering 
of intrinsic defects play the leading role. For instance 
paramagnetic samples can be turned ferromagnetic by 
prolonged exposure to Ti vapours, which produce high 
concentrations of Vcf^. 



FIG. 5: Proposed phase diagram for (Zn,Co)0 as a func- 
tion of the relative concentrations of Co and CoV. The yellow 
area is the blocked superparamagnetic (SPM) region, where 
both magnetic moment and hysteresis can be detected. The 
dashed line at small x delimits the region where most of the 
experiments are conducted. Finally the FM region is parti- 
tion into two regions depending on whether or not anomalous 
Hall effect can be detected. 

lows one the detection of both a net magnetic moment 
and hysteresis at room-temperature. For larger x a long- 
range FM ground state emerges, which however is lim- 
ited by the short range AFM interaction of the Co ions. 
Therefore we predict either a frustrated antiferromagnet 
or a spin-glass for > and <C We 

emphasize that bulk measurements (hysteresis or mag- 
netization) can hardly distinguish between the FM and 
the blocked superparamagnetic phase, and more local 



Methods 

Density functional theory calculations have been per- 
formed by using our newly developed approximate ver- 
sion of the self-interaction correction scheme (ASIC^f^. 
This is implemented numerically in the localized basis 
set code Siest^^. Here the ASIC scaling parameter 
a was set to 1/2, which is the value appropriate for 
describing the electronic structure of mid- to wide-gap 
semiconductors^^. In all the simulations we have consid- 
ered unit cells ranging from 128 to 256 atoms. The basis 
set used was as follows: Zn: DZ-s, DZ-p, SZ-d, O: DZ- 
5, DZ-p, SZ-d, Co: DZ-5, DZP-p, DZ-d (SZ=single zeta, 
DZ=double zeta, DZP=double zeta polarized^^). The 
grid cutoff (equivalent to plane wave-cutoff) was 650 Ry 
and we have considered 18 /c-points in the full Brillouin 
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zone for the 128 atom cell and appropriate scaling for 
other cells. Standard conjugate gradient geometrical re- 
laxation was performed until the forces were smaller than 
0.04 eV/A. 

Monte Carlo simulations were performed with the 
Metropolis algorithm as implemented in a home-made 
package. Typical simulation cells for the two center 
Heisenberg model contain between 1,600 to 5,000 mag- 



netic ions and we always use periodic boundary condi- 
tions. The systems are equilibrated until an initially 
AFM and FM replica have converged to the same value 
(typically 10,000 steps) and then the Monte Carlo mea- 
surements are taken by sampling 60,000 new steps. Dis- 
order averages are taken over 32 different samples at each 
concentration. Tq for the case of = 7% is evaluated 
from the specific heat (see supplementary materials). 
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